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We investigate the effect of modified gravity with screening mechanisms, such as the chameleon or 
symmetron models, upon the structure of main sequence stars. We find that unscreened stars can 
be significantly more luminous and ephemeral than their screened doppelgangers. By embedding 
these stars into dwarf galaxies, which can be unscreened for values of the parameters not yet ruled 
out observationally, we show that the cumulative effect of their increased luminosity can enhance 
the total galactic luminosity. We estimate this enhancement and find that it can be considerable 
given model parameters that are still under experimental scrutiny. By looking for systematic offsets 
between screened dwarf galaxies in clusters and unscreened galaxies in voids, these effects could form 
the basis of an independent observational test that can potentially lower the current experimental 
bounds on the model independent parameters of these theories by and order of magnitude or more. 
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I. INTRODUCTION 

Over the last decade or so, there has been a plethora of 
evidence that the universe is accelerating (see, for exam- 
ple [1] and references therein). It is well known that all 
observed matter that interacts via gravity as described 
by general relativity (GR) feels an attractive force. This 
should cause the universe to decelerate, or at least re- 
main stationary. Many of the attempts to explain this 
phenomena have tried either to introduce some exotic 
form of matter in the form of dark energy or to propose 
a modified theory of gravity. Many of these new forms of 
matter take the form of a scalar field and Weinberg [2 
showed that any modification to general relativity neces- 
sarily introduces a new degree of freedom and so these 
two different approaches are, on some level, equivalent. 
One would generically expect any new scalar degree of 
freedom to couple to ordinary matter, which invariably 
leads to the introduction of an additional, or fifth, force 
which alters the strength of the gravitational attraction. 
These forces are absent in general relativity, whose pre- 
dictions have been tested to high precision over the years 
in both the laboratory and the solar system (see [3] for 
some examples). In particular, these tests have placed 
very stringent constraints on any additional forces. Ei- 
ther they must be short ranged ( < 0.1 mm) or their 
strength must be weak compared to gravity. One is then 
naively led to conclude that these experiments rule out 
any theories of modified gravity that predict any non- 
negligible modifications to gravity. 

It has recently been realised, however, that these ex- 
periments have all been performed in the neighbourhood 
of our solar system and hence they do not rule out modi- 
fications to gravity where fifth forces are active over large 
distances but are suppressed on small scales so that these 
local constraints are satisfied. Theories such as this are 
said to possess screening mechanisms and regions where 
these mechanisms act to suppress the force are referred 
to as being screened. Such a mechanism was proposed 
by Khoury and Weltman in 2004 ;4] in the form of a 



scalar-tensor theory with a specific form of the scalar- 
gravity coupling. This chameleon mechanism, which is 
also present in F{R) theories [F, includes a coupling of 
the scalar to matter resulting in its properties becoming 
dependent on the ambient density. Ultimately, this has 
the effect of screening the force in high density regions of 
space-time (such as our solar system) and recovering the 
force predicted by general relativity. Since then, a num- 
ber of different mechanisms based on similar principles 
have emerged in the context of the symmetron [B] and the 
environmentally dependent dilaton [7]. There has been 
an intense amount of work developing these theories and 
in looking for observational tests. This has proved dif- 
ficult since the very mechanism that makes the theory 
viable and evades solar system constraints tends to sup- 
press any interesting observational features. Nonetheless, 
these theories can lead to subtle astrophysical effects such 
as modifications to galactic dynamics [SI [S] . 

In this paper, we present a new, model independent 
observational effect of modified gravity with screening 
mechanisms. Dwarf galaxies in voids can be unscreened 
in theories with values of the parameters not yet ruled out 
by observation whereas those in clusters are screened by 
the presence of their companions. It is possible then, that 
the stars within these galaxies are partially unscreened 
and so feel a stronger gravitational force in their outer 
layers. In order to support themselves against this in- 
creased gravitational force they must burn their fuel at a 
faster rate and hence radiate more energy in a given time 
period. As a result, their luminosity is greater and their 
lifetime is reduced. These dwarf galaxies are then more 
luminous than their screened counterparts in clusters due 
to cumulative effect of this increased luminosity in their 
constituent stars. We would then expect that one could 
construct a statistically based observational test of these 
theories by looking for systematic offsets between dwarf 
galaxies situated in voids and those in clusters. Recently, 
the effect of chameleon theories upon the structure and 
properties of red giant (RG) stars, which can be resolved 
in other galaxies, has been investigated numerically by 
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We analytically investigate this effect in order to gain 
physical insight into the the effects these theories can 
have on both stars and dwarf galaxies. For this reason, 
we consider only main sequence (MS) stars since they are 
far simpler in their structure and dynamical processes 
than their post MS counterparts, postponing a full nu- 
merical investigation for future work These theories 
can be described in a model independent manner using 
two parameters and there is currently some debate as 
to the observational constraints on these parameters de- 
pending on whether one demands that the galaxy is self 
screened or whether it is screened by other galaxies in the 
local group. Importantly, the effect presented here dif- 
fers from those used to place these current bounds and 
so could form the basis for an independent constraint on 
values of these parameters that may not yet be ruled out 
observationally. We have hence consider the entire range 
of parameters that are currently under debate. 

Main sequence stars are more readily screened than 
larger stars that have evolved off the MS and therefore 
an analysis including the effects on these stars could be 
used to probe parameters and order of magnitude (or 
even more) lower than have currently been tested and 
this is investigated in a followup work 

This paper is organised as follows. In section [ll] we 
present the model and derive the screening mechanism, 
elucidating this with common examples for the unfamil- 
iar reader in appendix A. In section [TiT] we calculate the 
effects of modified gravity on stellar structure and cal- 
culate the luminosity enhancement alluded to above and 
estimate the reduction in the main sequence life time. 
We also present some numerical calculations in order to 
reinforce these results although we shall not use them 
further in this work. We then go on to embed these stars 
in dwarf galaxies in section |IV| and estimate the galactic 
enhancement in unscreened galaxies due to this increased 
luminosity and shorter lifetime. We discuss our results in 
section [V| and speculate on other effects related to these 
phenomena that may provide additional tests of these 
theories. Throughout this work we use units such that 
h = c = 1 and use the metric signature (— , +, +, +). The 
reduced Planck mass is 



where is a scalar field with self-interactions given by the 
potential V{(j)). Here, is the matter action, are the 
matter fields, which are minimally coupled to the confor- 
mally rescaled metric g^i, = {(j)) g . The metric g^jy is 
known as the Einstein frame metric whilst g^i, is referred 
to as the Jordan frame metric[26]. Quantities (such as 
the covariant derivatives, energy-momentum tensor etc.) 
associated with each metric are distinguished by the use 
of tildes, with the indices raised and lowered using their 
respective metrics. The function A is known as the cou- 
pling function and is assumed to be a weak function of (j), 
i.e. A{(j)) 1 + 0{(l)/Mp\) so that metric perturbations 
in each frame are small. 

Matter fields couple minimally to the Jordan frame 
metric and hence it is the energy-momentum tensor in 



this frame, T^" 
served V^f.^;'' = 
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In the Einstein frame the energy-momentum tensor is 
given by 



2 SS„ 
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Since T^"^ is explicitly sourced by the scalar field, the 
energy-momentum tensor in the Einstein frame is not 
conserved 



d\n A{(j)) 



ray 
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where denotes the trace (with respect to g^^) of the 
energy-momentum tensor. The relation between the Ein- 
stein and Jordan frame energy-momentum tensor is 



(4) 



In this paper, we will work exclusively in the Einstein 
frame. In this frame there is an explicit coupling of the 
fifth force carrier (j) to the matter tensor in its equation 
of motion [27l 
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dKff 
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II. MODIFIED GRAVITY WITH SCREENING 
MECHANISMS 

A. The Model 



Meanwhile, varying the action eqn. ([T| with respect 
gives the modified Einstein equation 

= R". - IrS". = A/pi^ [T^, + T/,] , (6) 
where the scalar energy-momentum tensor is 



Theories of modified gravity where the extra degrees of 
freedom are "screened" on small scales may arise through 
the scalar-tensor action 



M, 
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^-Viq^) (1) 



^(V0)2 + !/(</)) 



(7) 



From this, it is clear that the scalar field (p experi- 
ences an effective potential Vcff (</>; — T^) which depends 
not only on but also the trace of the matter energy- 
momentum tensor T^. In this work, we will assume that 
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the matter is non-relativistic dust as in most astrophys- 
ical systems. The usual practice is to define the matter 
density via T„ = — p, however, as we have seen above, 
this quantity is not conserved in the Einstein frame. 
Now if matter is non-relativistic in the Jordan frame, 
p — —g^^T}!^ , then it can be shown that the quantity 
TJ^^ = A^{(I))T^^ is covariantly conserved in the Ein- 
stein frame provided that it is also non-relativistic in that 
frame [28 . If one then defines the energy density in the 
Jordan frame, p = —T„i, which is covariantly conserved 
in this frame, then one can showfT^ that the 'density' 
associated with Tm — ~A^{<f)p is (non-relativistically) 
conserved in the Einstein frame. For this reason, we shall 
work with this quantity as the density and not the one 
defined using T^. From here on we shall use the symbol 
p to denote this quantity and refer to it as the matter 
density. Since, to leading order, A{(j)) « 1 these two 
quantities do not differ from each other significantly, al- 
though one should bear in mind the discussion above. 
With this definition, the effective potential is 



VMp)^V{<t>) + pA{4>). 



(8) 



The geodesic equation for a test particle made up of 
'^i in the Einstein Frame then gives us the modified force 
law in the non-relativistic weak field and static limit [12] 



pi 



(9) 



where the (j) field explicitly appears as a source term as 
expected. The quantity 



d^(^) 



(10) 



then characterises the strength of this force and should 
only be weakly dependent on <j) if A{(j)) is not to differ 
too greatly from 1. 



B. The Screening Mechanism 

In this section, we review the physics behind the 
screening mechanism and derive an estimate for decid- 
ing whether an object is screened or not. 

The qualitative features of the screening mechanism 
are as follows: Suppose that the effective potential [8] pos- 
sesses a minimum. The density dependence of the poten- 
tial has the effect that this minimum occurs at different 
values of in regions of differing densities. If the field 
value in high densities is such that the additional force 
([9]) is negligible then any fifth forces will be screened. 
The form of eqn. Q reveals that there are two possible 
methods of accomplishing this. Firstly, the mass of small 
oscillations about the minimum is proportional to the 
density. If the density is high enough such that this mass 
is large compared with the inverse of the length scale as- 
sociated with the high density region then the field will 



only vary from its minimum by a small amount over the 
entire region and so the fieled gradient will be negligible. 
This can equivalently be thought of 8jS 8j Yukawa sup- 
pression. This is the mechanism employed by chameleon 
theories^. Secondly, if /3((/)) is close to zero in high 
density regions then the additional force will again be 
negligible. This is the mechanism employed by the sym- 
metron [13] and the environmentally dependent dilaton 
[7] .In principle, one could construct a theory where both 
of these effects are present but, to date, only theories 
where one of the effects are at play have been consid- 
ered. Of course none of these mechanisms can operate if 
the field cannot reach the value which minimises the ef- 
fective potential and so only regions where this is possible 
can effectively screen the fifth force. We thus require the 
field to minimise its effective potential throughout the in- 
terior of any high density region if the screening is to be 
effecient; this is illistrated in fig. (II B I. In appendix A 



we provide a brief overview of each of these mechanisms 
for the unfamiliar reader. 

With the above considerations in mind, we proceed to 
elucidate the screening mechanism in detail. We consider 
a homogeneous cosmological background with density po 
and assume that sufficient time has passed such that the 
scalar field has relaxed to the minimum of the effective 
potential 0o in such a background. The scalar in this 
case possess no dynamics, and hence carries no fifth force 
while its potential Vcs{(j)o) can act as a cosmological con- 
stant term which is small. We will use the subscript 
to define quantities in the cosmological background, i.e. 
/3o = /5(0o) and toq = (PV {(j)^) / d(jP from now on. 

Next, consider a small spherically symmetric overden- 
sity perturbation 5p{r) of radius R in the matter density 
around this background, which induces a perturbation in 
the field 

= 00+^0. (11) 

The field inside the overdensity then attempts to relax to 
its new minimum and the resulting field gradient gener- 
ates a fifth force on a test particle (see eqn. (jo])). If the 
perturbation is sufficiently large however, the field would 
reach its new minimum inside the perturbation resulting 
in a constant radial profile towards the centre and hence 
no extra fifth force is felt in this region. This central idea 
behind the screening mechanism: if the perturbation is 
large enough that the field can reach its new minimum 
(where, as discussed above, the fifth force is negligible) 
efficiently then there is only a small field gradient due 
to the interpolation between the two minima. If, on the 
other hand, the perturbation is small then the field will 
not differ significantly from its value in the background 
and the resulting field gradients will give rise to appre- 
ciable fifth forces. Let us investigate under which the 
perturbation is large enough to screen any fifth forces. 

The same over-density 5p sources a perturbation in the 
Newtonian potential, $n, which produces a gravitational 
force 



(12) 
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(a)Screened 




(b)Unscreened 

FIG. 1: The field profile inside a spherical overdensity embed- 
ded in a homogeneous background. The field minimises its 
effective potential at a distance Vs from the centre. Far from 
the overdensity, the field asemptotically aproaches the value 
which minimises the effective potential in the background, 00. 
The region r > rs is known as the unscreened region and here 
the fifth force is 0(1). In the region interior to rs the fifth 
force is negligible and this region is said to be screened. 



The Newtonian potential, in turn, obeys the Poisson 
equation 

V^ci>N ^ AnGSpir) + ... (13) 

where we have added the eUipses ... = 0{Mpi ^(V^)^) + 
■^V{4>) to remind the reader that, although in the Ein- 
stein frame the scalar field does source the Newtonian 
potential, we neglect them since the gradient terms are 
second order (recall that we have assume that (j) is small) 
and its potential is small by construction. |29) 

In the static limit, where the dynamical time scale of 
the perturbation is far greater than the time scale on 
which the field relaxes to its minimum (this will always 
be the case in this work), the equation of motion, eqn. 



(Is]), for (j) is = dVee/d(j). Inserting the perturbation 



eqn. ( 11 ) into this equation gives 



V2(0o 



d_V 

d(f) 



/3(<^o + 54>) 



Af, 



pi 



(po + M(14) 



There are now two limits depending on the amplitude of 
the perturbation. In the presence of a small perturbation, 
we intuitively expect the field perturbation to be small 



1^01 <C 0. We can then linearise Eq. (14) to give: 



(15) 



where we have dropped terms proportional to d/3/ dcj) as 
we have argued earlier that /3 is only weakly dependent 
on (j). Furthermore on scales longer than the Compton 
wavelength of the scalar R <C 1/mo, the effective mass 
term is negligible leaving us with 



pi 



5p{r). 



(16) 



In the other limit when the perturbation 5p is large, the 
field will be able to attain it's new minimum and so the 
perturbation is also large. In typical theories of interest 
the field value at the minimum in high density regions 
is much less than that low densities and so |(50| sa <j)Q. 
At distances deep enough in the interior of the pertur- 
bation such that the field reaches its new minimum the 
derivative of the effective potential vanishes and (see for 
example |^) 



dV 



so that 



pi 







p{r) 



(17) 



(18) 



i.e. the field gets trapped inside the new minimum at 
some .screening radius r < rg. The field profile be- 
comes constant if r is decreased further leading to no 
fifth forces. We say that the region inside a sphere of 
radius is screened. In general, rg can be zero {fully 
screened), < rg < R {partially screened) or undefined 
{unscreened). Unscreened objects will feel the full effect 
of the modified gravity, while fully screened objects are 
blind to its existence. 

Combining these two limits gives an approximation for 



valid for r <^ mn 



Mp{r)/Mpi 




rg < r m 
r < rg 







(19) 



By noting that the density perturbation Sp is related 
to the Newtonian potential via the Poisson equation eqn. 
( 13 1, this equation for S(j) can easily be integrated to find 



the field profile in terms of the Newtonian potential $n- 
For > 0, this is 
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5^{r) « -00 + 2/3oMpi 



(20) 



where H{x > 0) = 1 and H{x < 0) = is the Hcaviside 
function. Taking the hmit of the above as r/rs oo 
(and hence /3oA^pi$Ar(r) = 5(j){r) — >• 0) we can find an 
expression for the screening radius, which is defined im- 
phcitly by a function of the constants and the Newtonian 
potential, 



2/3oMpi 



^$N(rs)-r,$^(r,) >0. 



(21) 



In principle, we can use the Poisson equation eqn. ( 13 ) 



in conjunction with some density profile p to solve for r^. 
Notice however, that since > 0, rg > and $Ar < 0, 
it is clear that there exist no solutions (i.e. < 0) when 



|$w(r, ^0)1 < 



2/3oMp" 



(22) 



In this case, the object is fully unscreened. It is impor- 
tant to note that Xo is completely defined by the param- 
eters of the theory in question. It follows then that if 
the surface Newtonian potential > xoi the screen- 
ing radius > and the object will become gradually 
screened. This result gives us a rough guide on whether 
an object is fully screened or not, simply estimate the 
depth of the gravitational potential at the surface of this 
object (including the mass from the surroundings and its 
own) and see if eqn. ( 22 ) is satisfied. Note that if = i? 
then the fifth force is completely suppressed and — 
everywhere. 

Now that we have the complete field profile we can calcu- 
late the fifth force explicitly. For spherically symmetric 
bodies we have 



d$> 



GM{r 



(23) 



where M{r) is the mass enclosed within a shell of ra- 
dius r, and so differentiating eqn. (20) we find the total 



force per unit mass (gravity and modification) in the un- 
screened region is 



F — Fn + F^ 

^ d$N ^ /3o d(p 
~ dr Mpi dr 

GM{r) ^ 

1 + ao 1 



(24) 
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where q;(0) = 2/3^((/)), which defines and effective gravi- 
tational force 



F 



G{r)M[r) 



(25) 



The effect of modified gravity is then to produce an ef- 
fective change in the gravitational constant G — ?> G{r) = 
G(l -|- Q!off (f)) in the unscreened region, where 



aoff(r) = ao 1 



M(rg) 
M{r) 



(26) 



One can think of the quantity 1 — M{rs) /M{r) as a scalar 
charge, the existence of which is the source of equiva- 
lence principle violation in such theories. We can see 
that completely unscreened objects will feel a fifth force 
that is a factor ao that of the Newtonian force. We may 
then describe all theories with screening mechanisms in 
a model independent manner using two parameters: ag 
and Xo- cko specifies the strength of the fifth force whilexo 
(known as /^o in the context of F{R) models) controls 
the extent of the screening and can take on a range of 
different values in any one specific class of models, de- 
pending on both the position of the minimum and the 
coupling. Currently, there is some debate on the value 
of Xo that is consistent with local observations. If one 
demands that the solar system is screened by the poten- 
tial of the local group galaxies then the lowest bounds 
come from cluster abundance statistics |14) . which yield 
Xo ^ 0{1Q^^). If, on the other hand, one demands 
that the sun is self-screening then, as explained above, 
we must take xo ^ ^^{Rq) ^ 0(10"^) [11]. 



III. THE EFFECT OF MODIFIED GRAVITY 
UPON STELLAR STRUCTURE 

In this section we describe how these theories with 
screening mechanisms alter the structure of stars com- 
pared with standard gravity. Before examining the de- 
tails, it is instructive to first consider the effects in a 
heuristic manner. If the star is self-screening then we 
expect any effects to be negligible (or at least unobserv- 
able). If, on the other hand, the star is at least partially 
unscreened then in some part of its outer layers, the grav- 
ity will be stronger than the standard case. Stars support 
themselves against gravitational collapse by hydrogen fu- 
sion in their core to release energy and provide an out- 
ward pressure. We would then expect these reactions to 
proceed at a faster rate in order to combat the increased 
gravity. This will have two effects. Firstly, the star is 
more luminous and secondly, it will deplete its fuel sup- 
ply at a faster rate and will hence have a shorter lifetime. 

This effect was first noted by previous work by Chang 
and Hui [10] who found that solar mass M = Mq stars 
in the Red Giant Branch (RGB) phase of their life have a 
large increase in its stellar envelope reduces the depth of 
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its self gravity potential, resulting in its outer "mantle" 
layers becoming unscreened. In this section, we show 
that even in stars on the main sequence, a substantial 
part of the star's interior can be unscreened. Indeed, we 
will show that r^/RQ « 0.3 for in solar mass stars if 
Xo ~ 10^^. We emphasise that these are 0(1) effects, 
however, they are also highly degenerate with other stel- 
lar properties such as metallicity and other environmen- 
tal effects. In this paper, we will argue that while it is 
hard to disentangle these effects in individual star sys- 
tems, the gross effects of such modified gravity stars on 
galaxies may be observable. 



A. The Theory Stellar Structure in Modified 
Gravity 

The standard approach to stellar structure is to con- 
sider static, spherically symmetric stars of mass M and 
radius R. In order for the star to remain static, the out- 
ward force due to the pressure generated by the nuclear 
reactions must balance the gravitational force F{r) (per 
unit mass). This condition is expressed through the hy- 
drostatic equilibrium equation (HSE) 



dP 
dr 



-F(r)p(r) 



(27) 



where p(r) is the density at a distance r from the centre. 
The radius R is defined to be the value of r at which the 
pressure becomes zero, although this can be ambiguous 
and in more realistic models a definition involving the op- 
tical depth is more appropriate. The mass M{r) enclosed 
by a shell of radius r (so that M — m{R)) may be found 
by integrating the density over the volume enclosed by r, 
which leads to the continuity equation. 



dM 
dr 



47rr^p. 



(28) 



By considering photon propagation in the interior of the 
star, we obtain the radiative transfer equation, 



dr 3 K(r) p(r)L{r) 
dr 4a Anr^ 



(29) 



where L(r) and T{r) are the luminosity and temperature 
respectively at coordinate r. The quantity k is known as 
the opacity and represents the cross section per unit mass 
for radiation absorption; it is generally a function of the 
temperature and density. Finally, if energy is generated 
at a rate q{r) per unit volume by nuclear (or possibly 
other) processes then the luminosity gradient is deter- 
mined by the energy generation equation, 



dL 



= 4:nr^pq{r). 



(30) 



These are the equations of stellar structure, which de- 
scribe the relationships between the various thermody- 
namic and physical quantities which are required in or- 
der for the star to remain in a static configuration. As 



it stands, these equations do not close and so we require 
an equation of state for both the opacity and the pres- 
sure. We shall specify these later when we discuss scaling 
solutions and the Eddington standard model. 

Note that since Newton's constant, G, only appears in 
the HSE equation, the presence of modified gravity only 
changes that particular equation and not the others 150] . 
In other words, the total force acting on the star is the 
gravitational force and the fifth force given by eqn. (24 1. 
Eqn. (271 then requires that 



dP(r) 



GM{r) 



1 + ao 1 



M(r,) 
M(r) 



p(r). 



(31) 

To find the gravitational force, we integrate the Poisson 
equation eqn. ( 13 ) and use the continuity equation eqn. 



( 28 1 to obtain the solution to the gravitational potential 



as a function of mass M(r) 
d$N _ GM{r) 



dr 



M(r) =47r [ r'^dr'p{r'). (32) 
Jo 



Finally, in order to close the system of equations, we find 



rg from Eq. (21), which after integrating by parts and 



using the Poisson eqn. ( 13 ) yields an implicit solution 
for rs 



4nG 



R 



rp{r) dr = xo 



2/3oA/, 



pi 



(33) 



The implicity of the solution for Tg means that we have 
to iterate to find the complete solution. We will do this 
in the next section for a polytropic star. 



B. Scaling Relations 

In the next section we shall derive the structure of 
partially screened stars but before doing so we can ex- 
tract all of the essential physics by considering the sim- 
pler case of completely unscreened stars where the effect 
of modified gravity is to change the value of Newton's 
constant by a constant factor G G{1 + uq). In this 
case, we can use dimensional analysis to find the scal- 
ing of the stellar luminosity with the gravitational con- 
stant at fixed mass. By replacing each quantity Yi where 
Y = {AI{r),P{r),p{r),T{r),L{r)} by its characteristic 
value (say stellar mass or central density) in the equa- 
tions of stellar structure and replacing the derivatives 
by ratios (remembering to account for the correct sign 
depending on whether a quantity increases towards the 
centre or surface) we can obtain scaling relations between 
these variables. For example, the hydrostatic equilibrium 
equation eqn. (27 1 becomes 



R 



(34) 



where Pc a-nd Pc are the central pressure and density 
respectively. This process may be repeated for the other 
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equations to find the complete set of scaling relations. 
In order to complete this analysis we need to specify the 
form of the pressure. We consider two different sources 
of pressure, namely that due to radiation, 



and that due to the gas, which we take to be ideal: 



^ gas 



(36) 



where fimn is the mean molecular mass, which we take 
to be constant. Low mass stars are supported mainly 
by gas pressure and in this case the luminosity scales 
as L (X G^M^. More massive stars are hotter and are 
therefore supported predominantly by radiation pressure 
with a mass-luminosity relation L cx GM. Physically, 
stars supported by radiation pressure absorb radiation 
in their interior to prove the outward force needed to 
prevent gravitational collapse and less reaches the sur- 
face. If we then consider two stars of equal mass, one 
screened and the other completely unscreened we find a 
luminosity enhancement of 



Xo=0 



Low mass stars 
High mass stars 



(37) 



in unscreened stars. Clearly for typical values of ag, 
which is 1/3 for F{R) gravity [S], this increase is of order 
1 and so this effect could constitute a potential observa- 
tional test of these theories. 



The Eddington Standard Model and Partially 
Screened Stars 



Stars on the main sequence are in general supported 
by a combination of radiation and gas pressures P{r) ~ 



-frad + Pgas- The sct of cquations eqn. (27)-(30), combine 



with a set of equations of state (and modified gravity) 
represent a highly complicated set of differential equa- 
tions which are difficult to solve in general. A complete 
description of stellar structure will require a numerical 
method, which we will present in the section after this. 

In this section, we will instead solve these set of equa- 
tions in a particularly simple approximation known as 
the Eddington Standard model. While this model does 
not describe true stars, it is a fair approximation to main 
sequence stars. More importantly, it captures the essen- 
tial physics and we can gain a lot of intuition and under- 
standing of the physics underlying such systems under 
the influence of modified gravity. 

In the Eddington Standard model, we make the as- 
sumption that the radiation entropy per unit mass is con- 
stant throughout S'rad = iaT^ /3p — const. This greatly 
simplifies the equations by decoupling the HSE from both 
the radiative transfer eqn. ( |29| and the energy generation 
eqn. ( 30 ) . This assumption constrains T to be a function 



of p alone and therefore reduces the number of indepen- 
dent degrees of freedom in the system. Furthermore, by 
writing Prad = (1 ^ b)P, Pgas = bP where b measures 
the fraction of gas pressure support and assuming that 
Pgas is given by the ideal gas law, the equation of state 



(35) becomes polytropic i.e. P = Kp^l'^ |3T], as follows 



P = P 



rad 



K{h) 



3 / fc y {x-h) 



Kb 



1 1/3 



PT 



64 



(38) 
(39) 



Hence we simply need to solve the HSE to obtain P(r) 
and p{r) which would then fully describe the star. In ad- 
dition, the power law dependence of p means that the so- 
lution, in general without modified gravity, is self-similar 
(i.e. one can describe using dimensionless coordinates). 
The rescaled equations in this case is called the Lane- 
Emden equation. On the other hand, the presence of a 
screening radius weakly breaks this self-similarity, re- 
sulting in a modified form of Lane-Emden equation which 
we will solve below. Note that while b depends on the 
properties of the star (and modified gravity), it is a con- 
stant for a given star in the Eddington Standard mode, 
hence the space of Eddington Standard model solutions 
is spanned by this single family of parameters, b. 

Once a solution p{r) and P{r) is obtained (and hence 
r(r) by the assumption of constant entropy density), we 
substitute it back into the radiative transfer eqn. (291 to 



obtain the total luminosity at the star surface where r — 
R. Further assuming that the opacity is a constant and 
given by the electron scattering opacity taking Kop(P) ~ 
Kos we obtain the luminosity in terms of Ofsurf = ctoff (P) 
and M: 



L 



4^c(l-&)(l + asurf)GM 



where 



ao 



1 



M(r,) 
M 



acff(M). 



(40) 



(41) 



Before proceeding to solve the Lane-Emden equation, 
we pause to discuss the nature of these approximations. 
The decoupling of the energy generation and radiative 
transfer equations from the mechanical structure equa- 
tions have allowed us to close the equations. This means 
that we do not include any information about the nuclear 
processes occurring in the core, nor do we account for any 
chemical composition or gradient which may change over 
time. We have also used a constant opacity, which in 
general has a temperature and density dependence and 
have therefore ignored the effects of bound-bound and 
bound-free transitions, as well as H~ effects in the outer 
layers. 

Finally, we have not included the effects of convec- 
tion, which tends to be important in the outer layers of 
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stars. Nonetheless, the Lane-Emden solution is a good 
approximation to stars near the zero age main sequence 
(ZAMS). For this reason, we will not include the radial 
dependence of the fifth force in the following, but will 
model it by scaling Newton's constant by a constant fac- 
tor, G — > G(l + Q!o). When we come to the full numerical 
problem in later sections we will include the full radial 
dependence as well as a complete theory of chemical evo- 
lution and convection and so will be able to make far 
more realistic predictions for any evolutionary stage in 
the stars life. 

Let us now solve the equations. What follows is a tech- 
nical derivation, which the reader may not be interested 
in and is not required to understand the rest of the paper. 
Those who are not can safely skip to the next section. 

Since polytropic solutions are spanned by 6, our goal is 
to solve for b given ag and the stellar mass M. The Lane- 
Emden equation can be obtained by first defining a set of 
dimensionless variables - define P, p and T at r = to 
be Pc^ Pc and Tc respectively as the central values. Then 
we can define a dimensionless radial coordinate ^ by: 



Pr 



■kGpI 



1/2 



where Tc is the characteristic length-scale of the star and 
the dimensionless structure function 



(42) 



hence P = PcS^iO and p — pc9^{S). Inserting these 
new variables into the HSE eqn. ([27]), ignoring the radial 
dependent term M{r)/M(rs), and differentiating it then 
gives the modified Lane-Emden equation for ^(^): 



1 

e 



_d 



d9 



rs <r < R, 
r < rs, 



where the boundary conditions are the stellar radius R 
is defined by P{R) = and therefore 9{^r) — where 
— R/tq. Since we have scaled out the central pressure 
and density in terms of 9 we have 9{Q) = 1 and any sen- 
sible pressure profile should be smooth at the centre and 
so we require d0(O)/ d^ = 0. This second condition can 
also be viewed as a statement that there is no mass inside 
a shell of zero radius; setting M(0) = in eqn. ([27]) gives 
precisely this condition. With these two boundary con- 
ditions, we can proceed to solve the equation, matching 
up the solutions at r — rg. 

We can similarly rescale the screening radius and the 
total radius of the star by the characteristic scale i^s = 
Tg/rc and = R/r^; Ci? is determined by the require- 
ment that 9{£^fi) = i.e. p{R) = 0. The first derivative 
of the boundary condition is specified by 



09 



(44) 



«=€r 



Integrating this equation once (keeping rg and hence ^. 
fixed for the moment) and imposing continuity at re- 
quires us to further specify the matching condition 



de 



(45) 



Tg (and hence ^g) generally depend on both the mass and 
the radius of the star. It will turn out, however, that 
the Lane-Emden scaling constrains rg to depend on the 
combination xo/{GM/R). It then follows that to fully 
integrate the Lane-Emden equation we will need to find 
this quantity. 

Consider the case of standard GR {ao = 0), we have 
i^R — ~ 2.02. If on the other hand, the star is fully 
unscreened (ts = 0), then = and by rescaling ^ we 
have ujR ss 2.02(1 -I- ao)~^^'^. Between these two limits, 
we can define a function that interpolates between at G 
(0,ao) : 

/ _ X 2/3 

1 + + , (46) 

such that when rg = R, lur ~ lOs — Coi we have ab — 
0, while when rg = the star is fully unscreened and 
ctb = oiQ- The quantity ab is the analog of in Lane- 
Emden rescaled variables and is a function of b hence the 
its subscript. 

We can then find the mass of the star in terms of these 
variables: 

rR 



M 



in 



r'p{r)Ar = ATTp,4 ^9^00^, 



f 1 y/^ 


t^R ^ 




[ttGJ 


1 - 





(47) 



Inverting this gives 6 as a function of A/ through a mod- 
ified form of Eddington's quartic equation 



1 - b 



64 



= (1 



/ M 



cdd 



(48) 



with the Eddington mass 



1/2 



fcs y 



(49) 



which is defined in terms of fundamental constants and 
hence (Medd = 18.3/x~^M0) is independent of ao. Fi- 
nally, to close this set of equations we must specify ^g as 
a function of the modified gravity parameters (po and Pq. 
Using eqn. ( 19 1 we have 

/3o ^2 . 47rGQ!op(r) 



Mpi • " Mpi 

from which, using our knowledge that < 
we find that (ps — 0(7's) is given by 

Ms 47rGrgpc 



pi 



^R 

aoGM 
R 



1 + "0 



(50) 

as r — >■ oo. 



(r 



UJR + aoujs 
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The scaled screening radius, ^s, is then determined by 
4>s = -4>o and 



COR + aoWs 



R 



Xo 



2/3oAfpi GM GM/R' 



(51) 



This is the analog of the implicit solution for rg, eqn. 
(33), in Lane-Emden coordinates, and requires a similar 
iterative method to solve. If eqn. (51 1 has no solutions 
then Tg = = = and the fifth force is unsuppressed. 
This requires 



Xo 



GM/R 



> 1 - 



dlnf 



d6 



4.417, 



(52) 



for all ao- The numerical value 4.417 is a property of n = 
3 self-similar (i.e. completely screened) polytropes. Thus 
if GM/R > 0.226x0 the fifth force is at least partially 
suppressed. Finally, we can then find Ofsurf, which gives 
us the modification to the luminosity: 



1 



M{r, 
M 



ujR + ao 



(53) 



The procedure to find the luminosity increase is the 
following: we first solve the Lane-Emden equation nu- 
merically for a set of test matching radii . We then use 
this function 0{^',£,s) in conjunction with eqn. (|45| to it- 



eratively find using Eqn (471 and Eqn (51 1, and hence 
ab and Ofgurf ■ By repeating this for a range of different 
masses M and values of xo we numerically fit the lumi- 
nosity as a function of mass and xo- Since the screening 
radius is degenerate in xo and M/R, it is convenient to 
scan through the quantity X = xo/(4.417GM/i?) such 
that < X < 1. The relation between ah and ao can 
then be recast as ab{l + ab) = ao{l + ao)/(X; ao) and 
that ttsurf/ao = g{X]ao) where the fitting functions / 
and g take values between and 1. In doing this, we 
find f(X;ao) and g{X;ao) are independent of ao for 
ao ^ 0(1) and smaller and are respectively well fitted 
by the following functional forms 



^(1.94 + 0.79a; - 2.91x^ + l.lSa;^), 



9(.x) = \ - 




(54) 



In principle, one could vary ao to find the luminosity 
increase as a function of ao however in this work we con- 
sider only ao = 1/3, corresponding to F{R) gravity and 
the low background density limit of the environmentally 
dependent dilaton. The luminosity increase compared to 
the GR predictions is shown in fig[2j We can see that the 
enhancement is much greater in low mass stars in agree- 
ment with the simple scaling arguments presented above. 
Numerically, (1 + ao)^ = (4/3)^ « 3.16 and so we can 
see that for xo =10~^ and 10~^ the lowest mass stars are 
entirely unscreened. We also note that at 5xl0~^ there 



is an appreciable enhancement but at 10~^ we see that 
there is no significant effect. This is due to the fact that 
main sequence stars typically have Newtonian potential 
of order 10~^ and so, according the the arguments above, 
are almost entirely screened. After going off the main se- 
quence, the stars become red giants which have radii of 
10 to 100 times that of their main sequence phase and 
so even if they are screened in the main sequence, they 
can become unscreened in the red giant phase and show 
appreciable differences from GR [TO] . 
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»J Log[ 1 
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FIG. 2: The ratio of the stellar luminosity of a partially 
screened star compared with one that is completely screened 
as a function of stellar mass M. L^td is the luminosity of 
a completely screened (general relativity) star. From top to 
bottom: xo = (10"", 10"^ 5 x 10"^ 10"^). 



There is another important consequence of this in- 
creased luminosity. Although we have not included the 
effects of nuclear burning in our model, if the luminos- 
ity of a star in modified gravity is higher than the stan- 
dard case then the extra energy being radiated must come 
from an increased rate of nuclear fusion in the core. The 
main sequence lifetime is then shorter. The main se- 
quence lifetime can be well approximated by 



tms 



10 



M 
M^ 



L 



Gyr 



(55) 



and hence we can see that a fully unscreened star in F{R) 
gravity will go off the main sequence more than three 
times as quickly as one in general relativity. Ultimately, 
as we shall see below, this effect limits the luminosity en- 
hancement in unscreened galaxies since stars go off the 
main sequence before their enhanced luminosity can con- 
tribute significantly. This has important observational 
consequences since we then expect unscreened galaxies to 
have hosted many more generations of stars and therefore 
be more metal enriched than their screened counterparts. 
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D. Numerical Simulations of Modified Gravity 

The assumptions of the previous subsection are a 
good approximation to intermediate main sequence 
stars. However, as alhided to in the introduction to this 
section, making further progress and understanding the 
many degeneracies in stellar systems, requires accurate 
modelling of the stellar processes and knowledge of 
the complete evolution of the star from the pre-main 
sequence until its death. In practice, stellar evolution 
is determined using numerical evolution codes. In 
this section, we will take this approach for the case 
of modified gravity. As we have mentioned, modified 
gravity only changes the gravitational part of systems of 
equations, and hence we simply need to modify existing 
stellar evolution codes to make progress. We have chosen 
to implement the screening mechanism into the publicly 
available code MESA-star [16]. 

MESA is a one dimensional (in that it assumes spher- 
ical symmetry) stellar evolution code that numerically 
calculates the structure and evolution of stars by solving 
the stellar structure equations in conjunction with 
a fully consistent implementation of varying opacity 
laws, nuclear reaction networks, atmospheric models, 
convective theories, mass loss schemes and rotational 
dynamics. We refer the interested reader to the instru- 
ment paper [16j . We will describe our implementation 
of the screening mechanism into MESA and present 
some simple results in order to justify the analysis above. 

We start by noting that, given a density profile, the 
screening radius can be found using eqn. (331 given a 



specific value of xo- Using this result, we implement the 
modification of gravity into MESA using the following 
procedure. MESA divides the star into cells of unequal 
radius and assigns each cell a set of quantities such as 
radius, density, opacity etc. and we first integrate the 
quantity p from the surface of the star, cell by cell, until 
the condition l33l is satisfied. The radius of the cell where 
this is the case is marked as the screening radius. Next, 
we update the value of G in each cell using eqn. (26 1. 



The star is then allowed to evolve in time to the next 
model where this process is repeated. This process is 
accurate provided that the timescale over which the star 
evolves is not longer than the timescale on which the 
effective gravitational constant varies due to a change 
in the stellar structure. MESA automatically selects a 
timestep to reduce such numerical errors however it also 
has controls that allow the user to set the timestep and 
the number of cells and so an arbitrarily good precision 
can be achieved. 

In fig. ([3]) we show the Hertsprung-Russell (HR) di- 
agram for stars of one solar mass with initial metal- 
licity Z — 0.02 (solar metallicity) evolving from the 
ZAMS to the tip of the red giant branch for xo — 
5 X 10~^, 10~^, 10"'^ and the standard, unmodified case. 
It is clear from the tracks that stars that are less screened 



are indeed hotter and more luminous. We can also see 
that for Xo = 10~^ the tracks are identical along the 
main sequence but separate in the red-giant phase cor- 
responding to the Newtonian potential dropping due to 
the increased radius and the star becoming less screened. 
It appears that the red giant tracks all converge to a 
similar track. This is due to the fact that the Newto- 
nian potential of these stars is so shallow such that at 
the mantle they are unscreened to a very high degree, 
even at xo = 10"''. In fact, if one examines the tracks 
in detail then small differences can be discerned. Note 
that while the HR tracks for stars of similar masses may 
possess similar shapes, this does not mean that the other 
stellar properties are simply a translation in luminosity 
axis at the same age. 

To show this explicitly, we show the star's age and 
radius at three identical points along the stars main se- 
quence. Notice that it is clear that unscreened stars do 
indeed have shorter lives than their screened counter- 
parts. It also shows that the radii at the same evolu- 
tionary stage tend to be smaller as well. Physically, the 
extra pressure needed to support the star in modified 
gravity is produced by increased densities and tempera- 
tures over the standard case, demanding a more compact 
star. We have not shown the values for xo = 10^^ since 
these stars are almost entirely screened on the main se- 
quence and hence have nearly identical properties to the 
unmodified stars. Nor have we shown any information 
about the red giant phase stars. This is because the red 
giant phase is far shorter than the main sequence and so 
comparing quantities between stars that are screened to 
different extents is misleading. Despite the assumptions, 
it is clear that the luminosity increase described above is 
still present when more realistic models are used and the 
missing physics is accounted for. 



IV. GALACTIC PROBES OF MODIFIED 
GRAVITY WITH SCREENING MECHANISMS 

In this subsection we consider a fully unscreened 
galaxy composed of main sequence stars and calculate 
the luminosity compared with an identical unscreened 
galaxy. 

We have seen how the Newtonian potential determines 
the extent to which an object is screened. Ultimately, we 
would like to calculate the contribution from the stars to 
the galactic luminosity in the presence of modified gravity 
in unscreened galaxies. Treating an individual galaxy as 
a sphere of mass Mx and radius Rx, the characteristic 
potential is 



I* 



GM 



X 



(56) 



Now, dwarf galaxies typically have rotational velocities of 
the order 50-100 Kms~^ and hence, using the virial the- 
orem we have \GMx / Rx\ ~ 10~^-10~^. Dwarf galaxies 
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FIG. 3: The Hertsprung-Russell diagram for stars of one solar mass with initial metallicity Z = 0.02. The black line shows 
the tracks for stars in general relativity while the red, blue and green tracks correspond to stars in modified gravity with 
Xo = lO"'^, 10~® and 5 x 10~® respectively. The radius and age at the point where the central hydrogen mass fraction has fallen 
to 0.5, 0.1 and 10~^ is shown for each star except the xo = 10~^ case. 



are therefore an excellent probe of modified gravity theo- 
ries with Xo >10~^, at least one order of magnitude lower 
than current experimental bounds. Whereas these dwarf 
galaxies will not be self-screened, it is entirely possible 
that dwarf galaxies in clusters may be screened by the 
Newtonian potential of their neighbours. This will form 
the basis of any observational method of searching for 
these theories using stellar effects. 

The number of stars born with mass M within a galaxy 
is given by the initial mass function (IMF) $(M) = 
diV/ dM, i.e. the number of stars dN within mass range 
dM. In many stellar populations of, where we can re- 
solve individual stars, this relation is empirically found 
to be roughly universal |32| . For simplicity, in this work 
we use the Salpeter IMF [17], ${M) cx M'^-^s with 
0.08 < M < lOOAf ; leaving aside questions of whether 
this function is valid at very high and very low masses. 
We can then estimate the luminosity increase for an un- 



screened dwarf galaxy by integrating luminosity ( 40 1 us- 
ing the fitting functions (54) over the IMF. 



Before doing so however, we must account for the stars 
that have gone off the main sequence as the IMF only 
gives the number of stars born. We do this by making 
use of eqn. (55). If the age of the galaxy is Tago then 
we assume that all stars with tms > Tago contribute in 
their entirety to the luminosity whereas stars whose main 
sequence lifetimes are less than the age of the galaxy con- 
tribute a fraction TMs/^ago of their luminosity. We note 
that stars that have gone off the main sequence still have 
a luminosity enhancement in the red giant phase (and 
beyond), however we do not account for their contribu- 
tion here. This effect is accounted for in our analysis by 
including a factor /o(M;rage) where 



/o(A^;Tago) = 



1 

T MS 



Tms > Tago 

Tms < Tago 



(57) 
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so that the galactic luminosity is 

igal(Tago,Xo) = (58) 

.IOOM0 ^jy 

dM/o(M,Tage)L(M;xo)TT7- (59) 



We can then immediately see that stars whose main se- 
quence lifetime are shorter than the age of the galaxy do 
not contribute to the luminosity enhancement since the 
factor of L{M)^^ in the main sequence lifetime exactly 
cancels the factor in the integral. Normalising the in- 
tegrals so that the total luminous mass of screened and 
unscreened galaxies are identical (this is required in or- 
der to account for the fact that more stars have gone 
off the main sequence in the unscreened case), we have 
performed the integral ( 58 1 for dwarf galaxies of mass 



Mga,i =10^^Mq and age Tagc =13 Gyr for xo in the range 
10~^-10^^, a range which, as explained above, is not 
yet be ruled out observationally. Our results are sum- 
marised in table Table [l] reveals that the enhancement 



xo 


Luminosity Enhancement 


1x10-'' 


42% 


1x10"^ 


42% 


5x10-" 


29% 


IxlO-'' 


3% 



TABLE I: The luminosity enhancement in unscreened relative 
to screened dwarf galaxies as a function of xo- AH values were 
computed using qq = 1/3. 

is significant for xo ~ 10^^, which, as argued above, we 
would expect. The saturation around xo '^lO"^ is due 
to the effect of the decreased main sequence lifetime. If 
Xo is around this value then, as seen in fig. ([2]), the 
low mass stars are all completely unscreened and so in- 
creasing its value further cannot make them any more 
luminous. Hence, the luminosity enhancement saturates 
around this value. As noted earlier, red giant stars are 
still unscreened if xo ~ 10"^ and these stars tend to be 
far more luminous than those on the main sequence, see 
fig. ([3]) and so even at these low values of xo it is possible 
that there will be a significant effect, for example, galax- 
ies that have partially unscreened post main sequence 
stars and screened main sequence stars will produce dif- 
ferent spectra from those where all stars are partially 
screened. This may then allow observational bounds to 
be brought down by another order of magnitude or even 
further. 

It is clear that this analysis is only a first approxima- 
tion. We have not included the contribution from stars in 
any post main sequence evolutionary phase, nor have we 
accounted for any galactic effects such as outflows, dust 
obscuration and supernova feedback. In addition to this, 
the Eddington standard morel of main sequence stars do 
not account for important effects such as convection and 
metallicities. 



Despite these shortcomings, we argue that dwarf galax- 
ies located in intergalactic voids, which are not subject to 
the gravitational potential of nearby galaxies, should be 
unscreened whereas those in clusters are not and should 
therefore exhibit all of the effects described above. Of 
course one cannot hope to observe two identical galaxies, 
one screened, the other not, and simply measure the dif- 
ference in their luminosity. We would, however, expect 
there to be systematic differences between galaxies lo- 
cated in clusters and those in voids, of which the extent of 
screening can be estimated, which can be used as a test of 
these modified theories. By looking for these systematic 
differences independent constraints on xa in the range 
lO-^-lQ-® (or possibly lower) can be found. As men- 
tioned above, we also expect the unscreened galaxies to 
be more metal enriched, and so we expect there to be dif- 
ferences in other galactic properties such as the spectra, 
colour and magnitude. All of these possibilities are inves- 
tigated in a future work 11] using a complete numerical 
implementation of both stellar and galactic physics. 



V. CONCLUSION 

Scalar-tensor theories of modified gravity with screen- 
ing mechanisms are an appealing method of accounting 
for dark energy and the currently unexplained acceler- 
ation of the universe whilst still satisfying solar system 
tests of GR. These theories may be described by two 
model independent parameters, the enhancement factor 
of the additional force ao and xoj which parametrises 
how effectively objects can self-screen. In general, if the 
Newtonian potential is less than xo then it will be at least 
partially unscreened, ao is of 0(1) in any physically in- 
teresting theory and most experimental bounds on the 
theory come from looking at the values of xo that can be 
tolerated by experimental searches for fifth forces. 

If one demands that the Galaxy is self-screening then 
one must take a value ~ 10"^ but it is possible that our 
Galaxy is screened by the Local Group and in this case 
the bound is relaxed to 10"^. Dwarf galaxies have New- 
tonian potentials in the range IQ-'^ - 10-^ so that they 
are unscreened in all models satisfying these bounds; they 
are hence excellent probes of modified gravity theories. 
These rough estimates are based on the notion that main 
sequence stars are sufficiently dense to be self-screened. 
In the case where the galaxy is screened by the local 
group the bound of IQ-* does not come directly from this 
requirement but from cluster abundances and so tests us- 
ing astrophysical effects such as those presented in this 
work could improve this bound. 

In this work, we point out that main sequence stars 
are actually not fully self-screened; sufhciently low mass 
stars (and stars in the Giant branches) can and do feel 
the effect of modified gravity, provided that their envi- 
ronment is not screened. Hence, most stars in unscreened 
galaxies evolved very differently from those in screened 
galaxies. This means that we can use entire galaxies as 
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probes of modified gravity. This is a particularly crucial 
point - galaxies are many orders of magnitude brighter 
than stars and hence probe a much larger domain of space 
than stars, opening up the window for testing modified 
gravity to much greater length scales and allowing us to 
access available large data sets on galaxy spectra. 

Using a simple analytic model, we have calculated the 
effect of modified gravity with screening mechanisms on 
the structure of main sequence stars. There is a sig- 
nificant enhancement in the luminosity in models with 
Xo ^ 5 X 10~^. In addition to this, the star's lifetime 
can be severely reduced, with unscreened stars going off 
the main sequence more than three times as quickly than 
an identical screened star. An unscreened dwarf galaxy 
should then be far more luminous than one which is iden- 
tical but screened. While this is an idealised situation, 
there should be systematic offsets between the properties 
of unscreened dwarf galaxies in voids and those in clus- 
ters who are screened by the potential of other galaxies in 
close proximity. It may then be possible to use this offset 
as the basis for an observational test of these theories. 

We have estimated this enhancement by summing the 
luminosity enhancement of the constituent stars in the 
galaxy, appropriately weighted by the IMF and account- 
ing for stars that have gone off the main sequence since 
the galaxy was formed. In the most extreme cases, the 
total enhancement can be up to 40%. We have not inves- 
tigated the other effects of the reduced stellar lifetime. 
Stars that go off the main sequence faster inevitably die 
faster and so one would expect unscreened dwarf galaxies 
to have hosted many more generations of stars than iden- 
tical screened ones of the same age. We therefore expect 
systematic differences between the spectra and colour |33j 
in addition to the effects described here. These effects 
cannot be calculated in our analytic model and we inves- 
tigate them computationally in a followup work. 

In our analytic treatment, we have made several simpli- 
fying assumptions. Firstly, our stellar model does not in- 
clude the effects of metallicity or nuclear fusion. The stel- 
lar luminosity is degenerate with the metallicity and so 
a full analysis should account for stars of varying metal- 
licity within the galaxy. The lack of nuclear burning pre- 
vents us from calculating the dynamical evolution of the 
model and we therefore do not include the effects of dif- 



ferent reaction networks. Since the central temperature 
rises in order to combat the increased gravitational force, 
these networks may become significant in stars where 
they are usually negligible and this has a non-linear ef- 
fect on the stellar lifetime that we have not accounted for 
here. Other stellar processes such as convection may be 
important, however, they are too complicated to include 
in any analytic model. 

Secondly, our galactic model has not included non- 
stellar effects such as outflows and dust obscuration, 
which are likely to be important in determining the 
luminosity. Nor have we included the effect of post- 
main sequence stars, which are often more luminous 
than those on the MS. Since we expect these to be 
less screened than main sequence stars the effect of 
neglecting these is, in fact, to underestimate the galactic 
luminosity. In a follow up work, we shall numerically 
implement the screening mechanism into computa- 
tional models that fully account for the effects we have 
neglected here, some of which we have already presented. 

To check that our assumptions do reproduce the 
gross effects of modified gravity, we implemented the 
screening mechanism in a full numerical stellar evolu- 
tion code MESA. The results validated our assumptions 
and demonstrated the key features of modified gravity 
in stellar evolution : higher luminosity, shorter lifetimes, 
smaller radii, and higher effective temperatures. Using 
these numerical results of stellar evolution, we are con- 
structing realistic galaxies and their spectra with the goal 
of using them as precision probes of modified gravity. 



Acknowledgments 

First and foremost, we are incredibly grateful to Bill 
Paxton for tolerating and answering our many questions 
about MESA. We thank K. Koyama, A. Pipino, J. Chen, 
S. Hansen, D. Rudd, J. Zwart and especially L. Hui and 
P. Chang. ACD, EAL, JAS and DJS acknowledge the 
STFC, and EAL also acknowledges a NASA Astrophysics 
grant (09-ATP09-0049). 



[1] E. J. Copeland, M. Sami, and S. Tsujikawa, 

Int.J.Mod.Phys. D15, 1753 (2006), hep-th/0603057. 
[2] S. Weinberg, Phys.Rev. 138, B988 (1965). 
[3] C. Will, Pramana 63, 731 (2004). 

[4] J. Khoury and A. Weltman, Phys.Rev.Lett. 93, 171104 

(2004), astro-ph/0309300. 
[5] P. Brax, C. van de Bruck, A.-C. Davis, and D. J. Shaw, 

Phys.Rev. D78, 104021 (2008), 0806.3415. 
[6] K. Hinterbichler, J. Khoury, A. Levy, and A. Matas 

(2011), 1107.2112. 
[7] P. Brax, C. van de Bruck, A.-C. Davis, and D. Shaw, 



Phys. Rev. D82, 063519 (2010), 1005.3735. 
[8] F. C. Adams, JCAP 0808, 010 (2008), 0807.3697. 
[9] L. Hui, A. Nicolis, and C. Stubbs, Phys.Rev. D80, 
104002 (2009), 0905.2966. 
[10] P. Chang and L. Hui (2010), 1011.4107. 
[11] A.-C. Davis, R. Kotulla, E. A. Lim, and J. Sakstein 
(2011). 

[12] T. P. Waterhouse (2006), astro-ph/0611816. 

[13] K. Hinterbichler and J. Khoury, Phys. Rev. Lett. 104, 

231301 (2010), 1001.4525. 
[14] F. Schmidt, M. V. Lima, H. Oyaizu, and W. Hu, Phys. 



14 



Rev. D79, 083518 (2009), 0812.0545. 
[15] J. Khoury and A. Weltman, Phys. Rev. D69, 044026 

(2004), astro-ph/0309411. 
[16] B. Paxton et al., Astrophys. J. Suppl. 192, 3 (2011), 

1009.1622. 

[17] E. E. Salpeter, Astrophys. J. 121, 161 (1955). 
[18] P. Brax, C. Burrage, A.-C. Davis, D. Seery, and A. Welt- 
man, Phys. Lett. B699, 5 (2011), 1010.4536. 
[19] D. F. Mota and D. J. Shaw, Phys. Rev. Lett. 97, 151102 

(2006) , hep-ph/0606204. 

[20] D. F. Mota and D. J. Shaw, Phys. Rev. D75, 063501 

(2007) , hep-ph/0608078. 

[21] P. Brax, C. van de Bruck, D. F. Mota, N. J. Nunes, 
and H. A. Winther, Phys. Rev. D82, 083503 (2010), 
1006.2796. 

[22] R. Gannouji et al., Phys. Rev. D82, 124006 (2010), 
1010.3769. 

[23] D. F. Mota and H. A. Winther (2010), 1010.5650. 

[24] T. Damour and A. M. Polyakov, Nucl. Phys. B423, 532 
(1994), hep-th/9401069. 

[25] M. Gasperini, F. Piazza, and G. Veneziano, Phys. Rev. 
D65, 023508 (2002), gr-qc/0108016. 

[26] Our convention differs from that of [9] (HNS) by = 
5^^, A^{(j>) = fl{(p)-'^ and Mp(p = 4>. 

[27] The presence of the trace of the Einstein frame energy- 
momentum tensor Tm = gij.vT^ comes from a functional 
chain rule when we vary the action with respect with <j). It 
is clear since radiation is conformally coupled, T™"* = 
and hence do not enter into the equation of motion of (j). 

[28] Such "conservation laws" are not true in general, for ex- 
ample, if the matter is slightly relativistic or when ob- 
servers are highly boosted. 

[29] The innocuous nature of this equation tends to hide some 
important subtleties. Firstly, the Poisson equation in the 
Jordan frame possesses similar additional terms which 
can be neglected. In other words, the Poisson equation is 
approximately identical in both frames. This is a conse- 
quence of our imposition that A{cj)) ~ 1. Whilst in gen- 
eral this approximation does not have to be true for the 
screening mechanism to work, in practice its imposition 
allows us to treat non-relativistic matter fields in a sim- 
ilar way in both frames. In other words, dust behaves 
more or less like dust in both frames etc. although they 
are, of course, sourced differently by the energy density 
defined in their respective frames. Secondly, as an aside, 
note that second order terms with time derivatives do 
not appear in this equation due to our assumption that 
the system is static. 

[30] One can also imagine non-trivial couplings of the scalar 
to other degrees of freedom such as where F^^ 

is the electromagnetic field strength tensor, which would 
then change the other equations. We will not consider 
such couplings in this work although such couplings do 
indeed exist [18] . 

[31] This is a polytrope of index n = 3 which usually is a 
very good description of a radiation supported star. In 
addition, standard stellar structure notation uses /3 for 
fractional gas pressure support. Here, we use b instead in 
order not to confuse it with the modified gravity force 
coupling function. Other polytropic indices have been 
considered in the literature and our method and results 
generalise in a straight forward manner to these cases, 
however, for simplicity we present only the specialised 
case of n = 3 



[32] In principle, since the physics of gravitational collapse is 
expected to be encoded in the IMF, modified gravity may 
change its functional form. However, we do not attempt 
to investigate this here since the IMF is empirical and 
there is no known successful derivation either analytically 
or numerically of this function from known principles. 

[33] One could say that chameleons change the colour of 
galaxies. 



Appendix A: Specific Screening Mechanisms 

1. Chameleons 

First proposed by Khoury and Weltman in 2003 |15j , 
chameleon fields typically have a run-away potential such 
that F — >■ as (/) — >■ 00 (one standard paradigm being 
V{4)) = Af'*+"/0") and a monotonically increasing cou- 
pling function of the form 



A{<j>) 



(Al) 



(3 is usually assumed to be of ©(Afpi), corresponding to a 
gravitational strength force hovifever /3 3> 1 is not neces- 
sarily ruled out by current experimental bounds [THl HO] ■ 
In the standard case /3(0) = (34> where /? is a constant, 
however more complicated couphngs have been consid- 
ered [HHis]. 

With this choice of potential and coupling function, the 
field value that minimises the effective potential shifts to 
smaller field values in high density regions as shown in fig. 
( A 1 ). The chameleon screening mechanism relies on the 



effective mass being much larger in high density environ- 
ments compared with the cosmological value and hence 
the fifth force being incredibly short ranged. This can 
be seen to be the case in fig. ( A 1 ) . It is this "blending 
in with its environment" that gives rise to the particles 
name. In practice, this screening requires the chameleon 
to sit at the minimum of its effective potential through- 
out the majority of the interior of the over density and 
varying only in a thin shell near the surface. This is hence 
known as the thin shell effect. 



2. Symmetrons 

Symmetron fields [13] rely on local symmetry restora- 
tion as a method of screening and avoiding fifth force 
constraints. They are characterised by a Z2 invariant, 
symmetry breaking potential 



1/(0) 



with /i > and a coupling (also Z2 invariant) 

1 



A{ct>) 



1 



2M2 



(A2) 



(A3) 



which are both invariant under the symmetry transfor- 
mation — > — 0. We note that the potential and coupling 
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Veff 



cally), we have 








Veff 



(a) Small p 




Vi(f>) 



(b)Large p 

FIG. 4: The chameleon efTective potential (solid line) for small 
and large local densities. The dashed lines show the contri- 
bution from the potential and the dotted lines show the con- 
tribution from the coupling. 



'VEV 



±- 



'X 



(A5) 



The coupling function then gives the leading order 
coupling of small perturbations around this VEV to mat- 
ter as 



/3(0VEV) 



J pi'?' VEV 



(A6) 



Evidently, the coupling is proportional the local VEV and 
is hence zero in dense enough regions such that the sym- 
metry is restored but is strong when the density is small 
enough that the symmetry is broken. The symmetron 
parameters are chosen such that the symmetry is broken 
in our recent history in order to enter the epoch of dark 
energy domination at a similar time and avoid the coin- 
cidence problem. This requires fine-tuning of the mass 
scales M and jj. such that H'^Mpx' ~ ^^M^. In addition 
to this we require Af < 10~*Afpi in order to satisfy local 
gravity constraints. Unlike the chameleon and the en- 
vironmentally dependent dilaton (described below), the 
symmetron cannot account for the vacuum energy den- 
sity today. The vacuum energy at the symmetry broken 
minimum is of order HqM^ <C H^Mpi^ and so it is nec- 
essary to add a constant vacuum energy to the potential 
in order to reproduce the correct cosmic history. 



3. The Environmentally Dependent Dilaton 

The environmentally dependent dilaton [7^, is an en- 
vironmentally dependent generalisation of the Darmour- 
Polyakov mcchanism[21] that arises from string theory 
in the strong coupling limit. The starting point is the 
four dimensional low energy effective action in the string 
frame 



need only take this form in the vicinity of the minimum, 
and may have a more complicated general form at field 
values far away from this region. The potential eqn. ( A2 ) 



should hence be viewed as an effective field theory valid 
at low energies. The requirement that the symmetron 
force be comparable with that due to gravity in vacuum 
tunes the parameter A ~10~^^ and so the effective field 
theory is valid for all regimes of interest. The effective 
potential is then 



1 



(A4) 



and thus the symmetry is restored in sufficiently dense 
regions such that p > fi^M^. In high density regions of 
space-time the field's vacuum expectation value (VEV), 
'/'VEV = whilst in low density regions (i.e. cosmologi- 



S 



d*x 



A(0)-2- 



(A7) 



where (j) is the dilaton, Ig is the string length scale {Ig ^ 
with Ms the string scale), and we have allowed the 
coupling constants to contain some dilaton dependence 
not induced by the coupling function A. It has been 



argued [25 that the functions appearing in A7 take the 
following form in the strong coupling limit: 



bze-^ + 0{e~^^), 



2c| 



1 



(A8) 
(A9) 
(AlO) 
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with ci = Ms/Mpi > 1 and A - 0{l)-0{ci). Ap- 
plying this and making the conformal transformation 
5^11/ — {4') 9 fj-u we obtain the Einstein frame action 



S 



d^x^ 



where 



1 

A2 



dlnA 
d0 



(All) 



(A12) 



and V{4>) = A'^V{4>). Note that the field is dimension- 
less in this model, forcing us to change the definition of 
/3 relative to the standard case. The non-canonical ki- 
netic term modifies the equations of motion, which (if 
the contribution from the coupling constants is small) is 



1 



1 dKff 



2M, 



pi 



dip 



which then defines the effective potential 



(A13) 



(A14) 



We assume that the coupling function has a minimum at 
some (j) = so that near 0o we have 



(A15) 



with ^2 > 1 but A2{4> - 4>qY/2 < 1. These conditions 
are required in order to satisfy experimental constraints 
|7j. Minimising the effective potential and setting A{(j)) « 
1, we have 



/3(0min) =^2( 



p 



(A16) 



and thus in high density environments, the coupling func- 
tion tends to zero, whilst in low density environments 
/3 — >■ 1/4. Next, we must worry about the strength of the 
fifth force, which, using the non-canonical expression for 
the force enhancement is 



3A2/32(0)-|-r 



(A17) 



where we have used eqn. (A12). In this case, the theory 
is screened in high density environments, a(0niin) 
whilst in the opposite limit a{<j)) — >■ 1/3 for p — >■ where 
minimising the effective potential is equivalent to min- 
imising the potential alone. 



